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Abstract 

In this paper, we consider a leader-following consensus problem for networks of 
continuous-time integrator agents with a time- varying leader under measurement noises. 
We propose a neighbor-based state-estimation protocol for every agent to track the 
leader, and time-varying consensus gains are introduced to attenuate the noises. By 
combining the tools of stochastic analysis and algebraic graph theory, we study mean 
square convergence of this multi-agent system under directed fixed as well as switching 
interconnection topologies. Sufficient conditions are given for mean square consensus 
in both cases. Finally, a numerical example is given to illustrate our theoretical results. 

Keywords: consensus problems; multi-agent system; leader-following; stochastic 
system. 

1. Introduction 

In recent years, there has been an increasing research interest in the distributed co- 
ordination for multi-agent systems. This is partly due to its broad applications in many 
areas such as cooperative control of unmanned aerial vehicles, formation control [H 1261 [27] 
and swarming behaviors of social living beings [111 [T6| I30j . 

Consensus problems have a long history in computer science and formed the foundation 
of the field of distributed computing [12]. In consensus control, it is critical to design a 
decentralized network algorithm based on neighborhood information for agents to reach 
an agreement on their states, asymptotically in some sense. For a variety of consensus 
algorithms and convergence results we refer the reader to the comprehensive surveys \17\ 
[M] and references therein. Most researches in the previous literature assume the exchange 
of messages between agents is error-free. However, this is only an ideal approximation 
for real communication processes. Recently, consensus of dynamic networks with random 
measurement noises has attracted the attention of some researchers. In [8i|25j, the authors 
introduce time-varying consensus gains and design control schemes based on a Kalman 
filter structure. The decreasing consensus gain a{k) (where k is the discrete time instant) in 
the protocols is proposed in [5] to attenuate the measurement noises in a strongly connected 
circulant network. The analysis in [5] is generalized to strongly connected digraphs in 
[7] and digraphs containing a spanning tree in [6] by the same authors. The work in 
[10] deals with discrete-time average consensus problems in switching balanced digraphs 
under stochastic communication noises, while [9j investigates the continuous-time average 
consensus control with fixed topology and Gaussian communication noises. The authors in 
[13) treat a continuous-time leader-following consensus control under measurement noises 
with a constant state leader. 

In this paper, motivated by the above works, we consider a leader-following consensus 
problem for networks of continuous-time integrator agents with a time-varying leader in 
directed fixed and switching topologies. The control input of each agent is based on the 
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measurement of its neighbors' states and some estimated data of the leader which are both 
corrupted by white noises. We design a leader- following consensus protocol such that the 
leader has an underlying dynamics and some variables (e.g. velocity and acceleration) of 
the leader cannot be measured and every follower can obtain the measured information 
(e.g. position) of the leader only when they are connected with the leader directly. The 
collective behavior of self-organized groups of agents with active (or dynamical) leaders 
is one of the most interesting topics in distributed cooperative control. However, as [22] 
suggests, the extension of consensus algorithms from a constant reference to a time- varying 
one is non-trivial. Some related results can be found e.g. in [3l 1201 [3T]. where the systems 
considered are all deterministic and free of noise. 

Inspired by [5l|9l[l3], we introduce time- varying consensus gains in the followers control 
protocol to attenuate the measurement noises, which lead to a time- varying stochastic dif- 
ferential equation of the system. The state matrix of the equation is time-dependent and 
no longer a Laplacian matrix, and is neither symmetric nor diagonalizable. To implement 
the convergence study, we merge stochastic analysis and algebraic graph theory, by devel- 
oping a Lyapunov-based approach and addressing the Ito integral by the stopping time 
truncation method. Firstly, we derive a sufficient condition for the state of each follower 
to converge to that of the leader in mean square under fixed and directed interconnection 
topology. Then it is shown that the algorithm also render each follower track the leader 
in mean square under switching topology when the subgraph induced by the followers is 
balanced. 

The rest of the paper is organized as follows. In Section 2, we provide some preliminar- 
ies and present the leader-following consensus protocol. Section 3 contains the convergence 
analyses under directed fixed and switching interaction topologies. A numerical example 
is given in Section 4 and we conclude the paper in Section 5. 

2. Problem formulation 

Before we proceed, some basic concepts on graph theory (see e.g. [2]) are provided as 
below. 

Let G = (V, £, A) be a weighted digraph with the set of vertices V = {1, 2, • • • ,n} and 
the set of arcs £" C V x V. The vertex i in Q represents the ith agent, and a directed 
edge G £ means that agent j can directly receive information from agent i. The set 
of neighbors of vertex i is denoted by Mi = {j G V| G £}. A = (aij) G M"^" is 

called the weighted adjacency matrix of G with nonnegative elements and > if and 
only if j G Mi. The in-degree and out-degree of vertex i are defined as din{i) = Yl^=i 
and dout{i) = Z]j=i^ii) respectively. If din{i) = doutii) for i = 1,2, ••• ,n, then the 
digraph G is called balanced [18]. The Laplacian of G is defined as L = D — A, where 
D = diag((ij„(l), dini^), ■ ■ ■ , din{n)). A digraph G is called strongly connected if there is a 
directed path from i to j between any two distinct vertices i,j G V. There exists a directed 
path from vertex i to vertex j, then j is said to be reachable from i. For every vertex in 
digraph C/, if there is a path from vertex i to it, then we say i is globally reachable in G- 
This is much weaker than strong connectedness. 

Here, we consider a system consisting of n + 1 agents, in which an agent indexed by is 
assigned as the leader and the other agents indexed by 1, 2, • • • , n are referred as follower 
agents. The information interaction topology among n followers are described by the 
digraph G as defined above; and the whole system including n + 1 agents is conveniently 
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modeled by a weighted digraph Q = (V, f , A) with V = {0, 1, • • • , n} and 

/ ••• \ 

7 = 



aio ail • • • air 



)(n+l)x(n+l) 



where the lower right block submatrix of order n can be viewed as A. We define a diagonal 
matrix B = diag(6i, b2, - • • , bn) to be the leader adjacency matrix associated with Q, where 
h = ttio > and ftj > if and only if G Mi{G)- Here, Ni{G) is the set of neighbors of 
agent i in Q. 

The continuous-time dynamics of n followers is described as follows: 

Xiit) = Ui{t), i = 1,2, • • • ,n, (1) 

where Xi{t) £ M is the state of the ith. agent, and G M is the control input. The leader 
of this considered multi-agent system is described by a double integrator of the form: 

xo{t) = 9{t)vo{t) 
voit) = ao{t) (2) 
y{t) = xo{t) 

where g(t) : [0,00) (0, 00) is piecewise continuous, y{t) is the measured output and 
ao{t) is the input. We assume g{t) and ao{t) are known, that is, the dynamical behavior of 
the leader is precisely known (c.f. Remark 1). On the other hand, y{t) = X(){t) is the only 
data that may be gotten by the followers when they are connected to the leader directly. 
Since vo{t) cannot be measured, we have to estimate vo{t) in a distributed way during the 
evolution. The estimate of vo{t) by agent i is denoted by Vi{t), i = 1,2, ■ ■ ■ ,n. 

In our model, the ith agent receives information from its neighbors with measurement 
noises: 

yji{t) = Xj{t) + ajiUji^t), j G Mi, (3) 

yoi{t) = xo{t) + aoinoi{t), (4) 

where yji{t) {i G V,j G V) denotes the measurement of the jth agent's state Xj{t) by the 
ith agent. The {nji{t)\ j G V,i G V} are independent standard white noises and aji > 
is the noise intensity. 

A group of controls U = {ui \ i = 1, 2, • • • , n} is called a measurement-based distributed 
protocol [9], if Ui{t) G o'{xi{s),\Jj^j^ yji{s), < s < t), for i > 0, i = 1,2,--- ,n. 
Herein a{^\,\ G A) denotes the cr-algebra generated by a family of random variables 
{^A> A G A}. The so-called leader- following consensus problem is to design a measurement- 
based distributed protocol such that each follower's state will converge to the leader's in 
some sense as time goes on. 

Consequently, we propose the distributed control protocol which consists of two parts: 
• a neighbor-based feedback law: 

Ui{t) = h{t)(^ aij{yji{t) - Xi{t)) + bi{yoi{t) - Xi{t))'^ +g{t)vi{t), i = 1,2,--- ,n (5) 

where t > and h(t) : [0, 00) (0, 00) is a piecewise continuous function, called a time- 
varying consensus gain 
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• a dynamic neighbor-based system to estimate vo{t): 

Vi{t) = aoit) + jh{t) ( J2 ^ijiy'jiit) - ^iit)) + bi{yo,{t) - , i = 1, 2, • • • , n (6) 

where < 7 < 1 is some constant, and moreover y'jiit) and yoi(^) independent copies 
of yji{t) and yoi{t), respectively. In other words, we have 

y'jiit) = Xjit) + (Tjin'jiit), j £ Mi, (7) 

y'oM =Mt) + croin'oi{t), (8) 

where {n^j(t)| j E V,i G V} are independent standard white noises and independent with 
{n,i{t)\jeV,ieV}. 

The set of neighbors Mi of agent « in ([5]) and Q varies when the interconnection 
topology is switching and we defer the corresponding protocol formulation to Section 3.2. 
Remark 1. We take individual state Xi as scalar for simplicity in (OP and it can be 
extended to multi- dimensional scenarios as studied in 12(A \31f without much effort. For 
example, if Xi G M?, it can be thought as the position of agent i moving in a plane. 
Therefore, gvQ and gao + gvo are the velocity and acceleration of the leader respectively, 
which are known since the exact dynamics of the leader is assumed. 

Remark 2. We separate a factor g from the 'velocity term' of the leader in ^ in order 
to tone the decreasing consensus gain h, which appears to be a notable feature distinct 
from some kinds of uncertain environment (see e.g. \28l \29^). where a random term is 
directly appended behind the equation of the system. In such works, the consensus gains 
are supposed to have positive lower bound. 

Remark 3. From and it is clear that the designed protocol for the ith agent is 
indeed a measurement-based distributed protocol since it relies only on the state of itself 
and its neighbors. 

Let x{t) = - ,Xn{t))'^ and v{t) = {vi{t),--- Denote the ith row 

of the matrix A by Oi, and Sj := diag((Tij, • • • ,cr„j) for i = 1,2, ••• ,n. Hence S := 
diag(aiSi, • • • is an n x dimensional block diagonal matrix. Let no(t) = 

{noi{t),--- ,non(t))^ and ni{t) = (nij(t),-- - ,n„i(t))'^ for i = 1, 2, • • • , n. In addition, 
nQ{t) and ra-(t) can be defined in a similar way. The juxtaposed matrix Q := (B, S) is an 
n X n{n + 1) dimensional block matrix. Combining ([1]) with ([5]) and Q, we may write the 
protocol in a matrix form: 

= -h{t){L + B)x{t) + h{t)Blxo{t)+g{t)v{t) + h{t)QZ{t) 
= ao{t)l--fh{t){L + B)x{t) + -fh{t)Blxo{t)+-fh{t)QZ'{t) 

where Z(t) = (nj(t), (t), • • • ,n^(t))'^ and Z'{t) = {n'^^ {t),n\^ {t), ■ ■ ■ ,n'J{t)y are two 
n(n + 1) dimensional independent standard white noise sequences, and 1 = (I,... ,1)'^ G 
M". The system ([9]) may be further written in the form of the Ito stochastic differential 
equations: 

dx(t) = -h{t){L + B)x{t)dt + h{t)Blxo{t)dt + g{t)v{t)dt + h{t)GdWi{t) 
dv{t) = ao{t)ldt - -fh{t){L + B)x{t)dt + jh{t)Blxo{t)dt + -fh{t)GdW2{t) ^ ^ 

where Wiit) = {Wn{t),--- ,VFi„(t))^ and 1^2(0 = iW2i{t),--- ,W2„(t)f are two n 
dimensional standard Brownian motions which are independent with each other, and 

G ■■= diag(Jbl + ^j^^^a]^ajj , • • • , JW+T^^ZZ^U^j 
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3. Convergence analysis 

In this section we will give the convergence analysis of the system (|lUp and show 
that the state of every follower will track that of the leader in the sense of mean square 
convergence, that is, E\\x{t) — xo(t)l|| — > 0, as t ^ oo. Here || • || denotes Frobenius norm. 
Two different cases, fixed topology and switching topology, are explored. 
Remark 4. Mean square consensus protocols for stochastic systems are first introduced 
in 15] and then further studied by several researchers (e.g. l^\7^\^{Wi{Wil). Mean square 
convergence seems to be an important alternative for almost sure convergence in consensus 
problems under noisy environments. 

For a given symmetric matrix A, let Aniax(^) and Aniin(^) denote its maximum and 
minimum eigenvalue, respectively. To get the main result, we need the following assump- 
tions: 

(Al) The vertex is globally reachable in Q. 

(A2) There is a (5 > 0, such that ^ > ^(f^ + ^ for t > 0. 

Here, P is a positive definite matrix defined by Equation (jl3p . see below. 

(A3) /i(s)ds = oo. 

(A4) f^h\s)ds <oo. 

Remark 5. Assumption (Al) is imposed on the network topology, which is much weaker 
than strong connectedness. The technical Assumption (A2) roughly means that g is com- 
parable with the consensus gain h. Assumptions (A3) and (A4) are called convergence 
condition and robustness condition respectively in and which are standard assump- 
tions often used in the stochastic approximation fl5f . 

3.1. Fixed topology 

Let X* = X — xqI and v* = v — vqI. We then obtain an error dynamics of ([10]) as 
follows: 

de{t) = F{t)e{t)dt + G{t)dW{t), t>0 (11) 
where e{t) = {x* {t) , v* {t)f , W{t) = {Wi{t),W2{t)f and 

Here In denotes the n x n dimensional identity matrix. 
We will need a lemma for Laplacian matrix. 

Lemma l.[2l[23] The Laplacian matrix L of a digraph Q = (V, A) has at least one zero 
eigenvalue and all of the nonzero eigenvalues are in the open right half plane. Furthermore, 
L has exactly one zero eigenvalue if and only if there is a globally reachable vertex in Q. 
The main result in this section is given as follows: 

Theorem 1. For system ([I]) with the consensus protocols ([^ and if Assumptions 
(A1)-(A4) hold, then 

lim E\\e{t)f = 0. (12) 

t— »oo 

Proof. By Assumption (Al) and Lemma 1, we know L + B is a positive stable matrix, 
or in other words, —L — B is a stable matrix. From Lyapunov theorem, there is a unique 
positive definite matrix P E M"^" such that 

{L + BfP + P{L + B) = In. (13) 
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Let P = ^ ^ p ^ and define a Lyapunov function V{t) = (t)P£{t). Since 

< 7 < 1, P is a positive definite matrix. In fact, we have Amin(-P) = (1 — 7)Amin(-P) and 
Amax(-P) = (1 + 7)Amax(-P)- Utilizing Ito formula and pT|) . we have 

dV{t) = e^{t){PF{t) + F^{t)P)e{t)dt + tr(G^(t)PG(t))dt + 2e^ {t)PG{t)dW{t). 

Here tr(-) means the trace of a matrix. From the Lyapunov equation (jl3p . we get 

PFit) + F^it)P = -h{t) JJ, " Jili := -himt). 

Invoking the Haynsworth inertia additivity formula ;2l], Assumption (A2) and the positive 
definiteness of P, we know that Q(i) is positive definite with p := mint>o {Ainin(Q(i))} > 0. 
Thereby, we have 

AV(t) < -h{t) f^V{t)dt + /i2(t)tr(G^PG)dt + 2h{t)e^{t)PGdW{t). (14) 

Next we want to prove 

E [ h{s)e'^{s)PGdW{s) = 0, V < to < (15) 

For any given T > to > and K £ N, let Tj^''^ = T A inf{t > to| e^{t)Pe{t) > K}, which 
is a stopping time. By we have, for to < t <T, 



E{V{t At'^'')1^^^^.^o.t^) - EV{to) 

^* h{s)E{V{s A r^°'^)lj^^^^,Tj)ds + tr(G^PG) j' h\s)ds 



< P- 



<tr(G^PG) / h^{s)ds. 

Jto 

This implies that there is a constant Cto,T > such that 



E{V{t A r^°''^)l^^<^^o.Ti) < a„,T, V to < i < T. 
*o,T 



Since limx^oo t At'^''^ = t a.s., for to < i < by Fatou lemma, we derive 



sup EV{t) < Cto,T. 
to<t<T 

Accordingly, 

rt t-T 



e[ h'^{s)V{s)ds < sup EV{s)- [ h^{s)ds<(yo, V < to < t. 
Jto to<s<t Jo 



Combining this with 

ft ft 



eI h^{s)\\e^{s)PGfds<\\P\\\\GfE I h^{s)V{s)ds 

Jto Jto 
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yields (US} (c.f. [H]). 

Now employing (jl4p and (jlSp . we obtain 



- ^y(O) < /" /i(s)-Ey(s)ds + tr(G^PG) /" /i2(s)ds, V t > 0. 

Amax(-P) Jo Jo 

Thus, from the comparison principle |14| . 

EV{t) < EV{0)exp(^ / /i(s)ds) 

(P) ^ 

+ tT{G^PG) [ h^{s)exp( f h{u)du)ds. (16) 

Clearly, by Assumption (A3) the first term on the right hand side of (jl6p tends to 0, 
as t ^ oo. For any r/ > 0, by Assumption (A4), there exists some sq > such that 
h'^{s)ds < 7]. Hence, 

/i^(s)exp( ^ ^ / /i(n)dn)ds 

^ A^ax(P) Js ^ ' ^ 



= /\2(s)expf / /i(n)dn')ds+ / /i2(s)expf / hiu)du]ds 

<exp( [ h{u)du) [°h^{s)ds+ [ h^is)ds 

^ Amax(-P) J so ' Jo J so 

< exp ( ^ ^ / /i(n)dn j / /i^(s)ds + 77, V t > sq 

^ Amax(-P) J so ' Jo 

By Assumptions (A3), (A4) and the arbitrariness of r/, the last expression tends to zero, 
as t ^ 00. Therefore, (fTHj) yields limf^oo^^(i) = 0. Note that 

2< 



Amin(P) 

which concludes the proof. □ 

Remark 6. As is known, the solution to Lyapunov matrix equation may be obtained by 
using Kronecker product. Thus the positive definite matrix P involved in Assumption (A2) 
can be given explicitly. 

Remark 7. Theorem 1 implies that in the fixed topology, under Assumptions (A1)-(A4), 
the designed protocol guarantees the state of each follower tracks that of the leader in mean 
square. 

3.2. Switching topology 

In this section we deal with the convergence of the protocol under switching topology. 

Let a{t) : [0,oo) Sj-i = {1,2, ••• ,N} be a switching signal that determines the 
communication topology. The set is a set of digraphs with the common vertex set V 
and can be denoted as 7i = {Gi,G2, • • • ,Gn}, where is the total number of digraphs 
in H. Naturally, let Ga{t) be the subgraph of Q(j(t) induced by V. Thereby, we rewrite the 
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consensus protocols ([5]) and ([6]) as: 

Ui{t) = h{t)(^ a^j{a{t)){yji{t)-Xi{t))+bi{a{t)){yoi{t)-Xi{t))^ 

+ g{t)vi{t), i = 1,2,-- - ,n (17) 

and 

Vi{t) = ao{t) + ^h{t) 

•( «*j(^W)(yiiW-^*(*)) + ^^('^W)(yo*W-^*W))' i = l,2,---,n (18) 

jeM(<T{t)) 

where, Mi{a{t)) is the set of neighbors of agent i in the digraph Qa(t) formed by n fol- 
lowers; aij{a{t)) is the (i,j)-th element of the adjacency matrix of l^t -Bo-(j) := 
diag(6i((j(i)), 62(<7(i)), • • • , 6n(o'(t))) represent the leader adjacency matrix associated with 
Ga{t) such that bi{a{t)) > if and only if G ■N'iiQ a{t)) ■ 

In parallel with Section 2, substituting the protocols p!7|) . to the system ([T|), we 
can describe the system in the form of the Ito differential equations: 

r dx{t) = -h{t){L,(^ + B^^t))x{t)dt + h{t)B^^t)lxoit)dt + git)v{t)dt + h{t)G^^t)dWi(t) 

\dv{t) = ao(t)ldt-7/i(t)(ia{t) +5a(t))a;(t)dt + 7/i(t)B,(t)lxo(t)dt + 7/i(t)G'<,(t)dTy2(i) 

(19) 

where L^(^^^ is the Laplacian matrix of Ga{t) ) ^^'^ 

G.(,) :=diag( te(a(t))+ Y ^]AM^)),--- , te(a(t))+ Y 

V i6M(a(t)) V ieA/'n(a(t)) 

Let rr* = X — xqI and v* = v — vqI as in Section 3.1. We get an error dynamics of (jl9p 
as follows: 

de{t) = F^{t)e{t)dt + G^{t)dW{t), t > (20) 
where e{t) = {x* (t) , v* (t))'^ , W{t) = {Wi{t),W2{t))^ and 

- ( -"X'lfe;;, + b:;;,', ^*'„"" ) ■ - ^<') ( .g-::, ) - ^('x^-.'.- 

In the sequel, we show that under switching topology, the consensus protocols ()17p and 
p8|) ensure that each follower tracks the leader in mean square. We will use the following 
lemma. 

Lemma 2.^ Given t >0 and suppose the digraph Ga{t) ^-^ balanced. Then Lfj(i-^+B„^i-^ + 
{La{t) + -Scr(t))"^ positive definite if and only if the vertex is globally reachable in Ga{t) ■ 

The matrix i(j{t) + -^o-(t) + {La(t) + ^o-(t))^ plays a key role in the convergence analysis 
below. Define ^ := mint>o {Amm(^ff(t) + ^^(t) + (-^^o-(t) + ^a(t))^)}- 

Prior to establishing the main result, we present a condition analogous with Assump- 
tion (A2) in Section 3.1: 

(A5) There is a (5 > 0, such that > 2^(il^ij)^ + 5 for t > 0. 

Remark 8. It is easily shown that fi > under the assumptions of Theorem 2 below by 
exploiting Lemma 2 and the fact that Tt is a finite set. This validates the expression in 
Assumption (A5). 
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Theorem 2. For system ^^ with the consensus protocols |77| j and jTgjj . if for any t >0, 

Qa{t) is balanced, and vertex is globally reachable in Q^(t-), then under Assumptions (A3)- 
(A5), we have 

lim = 0. (21) 



Proof. Let / := f J"^" | . Obviously, we have Amin(^) = 1 — 7 and Amax(-^) = 

1 + 7. Hence / is a positive definite matrix by recalling < 7 < 1. Define a Lyapunov 
function V{t) = e^{t)Ie{t). 

By Ito formula and (j20p . we have 



dV{t) = e^{t){lF,{t) + Fj{t)l)e{t)dt + ii{Gl{t)IG„{t))dt + 2e'^ {t)IG„{t)m{t). 
Straightforward calculation yields 



IF,it) + Fjit)I 




(1 - 72)(L,(,) + + (L,(,) + i?.(,))^) -/„f| 
_r !?{*) o^r 



By using the Haynsworth inertia additivity formula [21j . Assumption (A5), we get that 
Qa{t) is positive definite with u := mint>o {Amin(Qo-(i))} > 0- 
Therefore, we have 



dV{t) < -h{t)-^V{t)dt + /i2(t)tr(G^(,)/G,(t))dt + 2h{t)e^{t)IG,^t)dW{t). (22) 



The remaining proofs are similar with those in Theorem 1 by noting that 7^ is a finite 
set, and hence omitted. □ 

Remark 9. From Theorem 2 we see that the designed protocol may guarantee the state 
of each follower tracks that of the leader in mean square even under the switching topology. 

4. Numerical example 

In this section, we provide a numerical simulation to illustrate the theoretical results. 

We consider a network consisting of four agents including one leader labeled by vertex 0, 
as shown in Fig. 1. The digraph in this figure is assumed to have 0—1 weights. With simple 
calculation, it is not hard to solve P from Equation ()13p and obtain Amax 

(P) = 0.9447. We 

take (Tjj = 0.1 for all i,j, h{t) = g{t) = -^(p^ and 7 = 0.5. Therefore, Assumptions 
(A1)-(A4) in Theorem 1 hold. 

The simulation results for the consensus errors x* and v* are shown in Fig. 2 and Fig. 
3 respectively, with initial value e(0) = (—2, 1.5, 3, 2, —1.5, —1)^. From Fig. 2 and Fig. 3, 
we can see that three followers can eventually follow the leader. 

5. Conclusion 

This paper studies a leader-following coordination problem for multi-agent systems 
with a time- varying leader under measurement noises. Although the state of the leader 
keeps changing and the measured information by each follower is corrupted by white noises, 
we propose a neighborhood-based protocol for each agent to follow the leader. We present 
sufficient conditions for each follower to track the leader in mean square under directed 
fixed topologies. Sufficient conditions are also provided when the interaction topology 
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is switching and the subgraph formed by the followers is balanced. Finally, numerical 
simulations are presented to illustrate the theoretical results. Topics worth investigating 
in the future include time-delay cases and the design of almost sure consensus protocols. 
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Figure captions 



Fig. 1 Directed network Q of four agents involving one leader. Q has — 1 weights. 
Fig. 2 Consensus error x* for the agents. 
Fig. 3 Consensus error v* for the agents. 
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